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Abstract 

We calculate the equation of state of the Gross-Neveu model in 2+1 dimen- 
sions at order 1/A^, where is the number of fermion species. We make 
use of a general formula valid for four- fermion theories, previously applied to 
the model in 1+1 dimensions. We consider both the discrete and continu- 
ous symmetry versions of the model. We show that the pion-like excitations 
give the dominant contribution at low temperatures. The range of validity 
for such pion dominance is analyzed. The complete analysis from low to high 
temperatures also shows that in the critical region the role of composite states 
is relevant, even for quite large N, and that the free-component behaviour at 
high T starts at about twice the mean field critical temperature. 
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Typeset using REVT^ 



I. INTRODUCTION 



To quantify the contributions of fundamental and composite degrees of freedom to the 
thermodynamics of strongly interacting matter is a very difficult task. 

On general grounds, one expects that at very low temperatures the light and relatively 
weakly interacting pions should dominate the partition function, whereas at high tempera- 
tures the relevant degrees of freedom are believed to be almost free quarks and gluons. 

Unfortunately, when studying chiral symmetry breaking and restoration in QCD, where 
the non perturbative regime plays a major role, we are not able, at present, to describe in 
a quantitative way, starting from the fundamental lagrangian, what happens in going from 
low to high temperatures. 

Quantitative results concerning the low-temperature behaviour of the fermion condensate 
and of the pion decay constant for massless quarks have been presented by Gasser and 
Leutwyler who make use of an effective lagrangian for pions, assumed to dominate the 
thermodynamics in this regime. 

Other studies, complementary to Ref. [|l| and to lattice results, have been generally per- 
formed by means of models which encompass some properties of QCD, with the nice feature 
of starting directly from the fundamental degrees of freedom. Nevertheless, they generally 
do not go beyond a mean field description, neglecting the role of bosonic fiuctuations. This 
approximation is clearly unsatisfactory for low temperatures, and also it only gives a very 
rough description of the transition. It is important to go beyond, at least to understand 
what "low" and "high" temperature mean, i. e. the regions where in QCD only pions or only 
quarks and gluons would respectively dominate, and also to give an improved description of 
the region in between. 

Recently, a certain number of studies have been devoted to this program, by means of 
large- expansions 0, N being the number of fermion species, carried out directly at the 
level of the functional integral describing the partition function. This procedure allows to 
preserve the desired symmetries from the beginning 0. 
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At finite temperature, this technique has been applied to the Nambu-Jona-Lasinio model 
and to the 1+1 dimensional Gross-Neveu model [§,0. In spite of the approximations 
involved (and mainly the fact that none of them starts from the fundamental QCD la- 
grangian), these studies are useful towards a better description of chiral symmetry breaking 
and restoration. 

Here we present a calculation in the Gross-Neveu model in 2+1 dimensions. The model 
exhibits chiral symmetry breaking for couplings stronger than a critical value, and it is 
renormalizable in the 1/A^ expansion [pHl^. The finite temperature extension of the model, 
in the mean field approximation, has been studied for instance in Ref. |l^,|T^. At order 1/A^, 



the zero temperature effective potential has been calculated in Ref. |15|, whereas the finite 



temperature gap equation has been discussed, by using the Ward-Takahashi identities, in 



Ref. \16 



We have calculated the finite temperature propagator of composite states which enters a 
general formula valid for four-fermion theories to calculate the finite temperature equation 
of state. This procedure had already been followed in studying the Gross-Neveu model in 
1+1 dimensions Here we perform, for Gross-Neveu in 2+1 dimensions, a full study of 
the equation of state and we quantify the role of the various contributions from low to high 
temperatures. 

We have studied both the version of the model with discrete chiral symmetry, which is 
spontaneously broken by mass generation, and the version with continuous chiral symmetry 
where, due to the effective 2-dimensionality at finite temperature, there is no true phase 
transition. Thus we have included a bare mass which avoids I. R. divergencies in this case. 

It is worth noticing that the full pressure at order 1/A^ is real for any temperature. This 
apparently trivial feature is not so obvious, since in general the zero temperature effective 
potential is not always real for small values of the order parameter, thus indicating either 
some instability or ill-definiteness of the theory at this level [^,|l^,|r^. The fact that the value 



of the effective potential at the absolute minimum turns out to be real for any temperature 
allows for a safe A^-expansion of the model. 



In conclusion, by exploring the approximation to the equation of state in the whole 
range of temperatures it comes out that the naive expectation that pions dominate the 
partition function at low temperatures is respected, as well as the free-fermion dominance 
far above the critical temperature. Nevertheless the results show that the pion dominance 
is limited for temperature ranging from zero to about 0.1 T^, for a reasonable value of the 
chiral symmetry breaking term. 

Furthermore, in the critical region, where nonperturbative effects are crucial, the mean 
field pressure is strongly modified by the correction, even for large A^. 

The rest of the paper is organized as follows: in Sect. |I| we recall the results for the 



renormalized zero temperature effective potential, at order In Sect. |T| we give the 

expression for the pressure at order according to Ref . and briefly recall the solution 
of the mean field gap equation. In Sect. |^ we comment the results concerning the discrete 
symmetry (Subsect. A) and the continuous symmetry (Subsect. B) versions of the model. 
Concluding remarks are in Sect. |V|, whereas in Appendix ^ we give the expressions for the 
inverse bosonic propagators entering in the total pressure. 



II. EFFECTIVE POTENTIAL AT ORDER 1/N 



In this paper we consider the four fermion model, in D = 2 + 1 dimensions, defined by 



the following Lagrangian |[l8| , p!9 | 



A 



A, 



2N 



2N 



(2.1) 



where ip is a multiplet of A^ four-components fermion fields, and 7^ is a traceless matrix 



which anticommutes with all 7^^ (/i = 0, 1,2) Notice that in dimensions D = 2 + 1, 

in order to allow for chiral symmetry, we cannot consider the two- components spinorial 
representation of the Lorentz group S0{2, 1), since there is no 2 x 2 matrix anticommuting 
with the 7^^. On the contrary, in the case of four-components fermion fields, there are two 
matrices, 7^ and 7^, which anticommute with 7°, 7^ and 7^ (see e. g. Ref. PD[ for details). 
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In the Lagrangian we have introduced two couphngs, A and A5, for the scalar and 
pseudoscalar part of the interaction respectively. In the following we consider the cases 
A5 = A and A5 = 0. Since the expressions for the effective potential and for the pressure for 
A5 = A can be easily reduced to those of the A5 = case by simply omitting the pseudoscalar 
contributions, in this section we derive all the formulas in the former case. 

For vanishing bare fermion masses, M = 0, the Lagrangian (|2.1| ) is invariant under the 
discrete chiral transformation 

(2.2) 

for A5 = 0, and under the continuous chiral transformation 

^ ^ e'^^'^ (2.3) 



for A5 = A. In the following, when considering the finite temperature extension of the 
model with A5 = A, the mass term M will be always retained finite, in order to avoid I. R. 
divergencies, which would otherwise occur, according to well known no-go theorems P,p!4|. 
Furthermore, since we are interested in looking for hints for QCD, where the quarks are 
massive, the exclusion of the case M = is not a severe limitation. 

The expression for the effective potential at order is P,P,p!2[] (a = M/A) 

where iT^Q j is the zero temperature inverse bosonic propagator given in Appendix ^ The 
previous expression is renormalizable, since four fermion theories are renormalizable in the 
expansion |]8|-[12|. At the leading order the U. V. divergent terms are proportional to 



ip"^, while at the next-to- leading order there are divergencies proportional to (f"^ and To 
renormalize the expression ( |2.4| ) we introduce a three-momenta U. V. cutoff A and write 

V(y^) = V"-'--(¥p) + ^^'5Z2{A) + \^\'5ZsiA) (2.5) 



To fix the finite part of the counterterms we follow the prescription of Ref. |T^, from which 
one gets 



j=(T,7r 



o7i Zn 



(2.6) 



(2.7) 



y 
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m-iTT'O + T; tan I — - 

2 y V2|<^| 
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(2.8) 
(2.9) 



where rriQ = 2 — tt/A is the minimum of Vo{(f) at a = 0. From this relation and from the 
study of the renormalization group /3-function it follows that the leading order of the theory 
has a fixed point X* = tt/2 (see Ref. fl^) which separates the two phases in which the chiral 
symmetry is broken (A > A*) and unbroken (A < A*) ||9,10,12]. The order 1/N leads to a 
correction of the critical coupling of the form 



A* 



A* 



2 



(2.10) 



valid for the zero temperature effective potential ( p.6|) with only the scalar term ( p.8| ) and 
also when the pseudoscalar ( |2.9|) is included . 



It is worth noticing that in studying the mean field approximation, it is not necessary to 
specify the value of tjiq or equivalently of the coupling A (supposed to be > A*). It is in fact 
simple to verify that the effective potential can be rescaled by mo by a re-definition of the 
fields and of the mass parameter a. This is no longer true at the order Thus, when 

studying the equation of state at order although the gap equation needs to be solved 
only at the mean field level, attention must be payed to choose a value for the coupling 



A > A*, namely consistent with chiral symmetry breaking at zero temperature for a = 0. 
When presenting our resuhs, it will be supposed that A = tt = 2 A* (and thus mo = 1) 
which is greater than A* for N > 2. 

A potential difficulty of the expressions in Eqs. ( |2.8| ) and ( |2.9| ) is that they are not always 
real for small values of ip (see also Anyway they are real around the mean field solution. 
This feature allows to evaluate the pressure in a scheme. 

Notice that all the variables in the previous expressions are dimensionless. They have 
been rescaled by the (arbitrary) renormalization scale which serves as a subtraction point 



m 



III. PRESSURE AT ORDER 1/N 

According to the expansion derived in Ref. ||^ the expression for the total pressure 
at order can be written as 

p(T) =Po(^o(T),T) + i y: nmnr) (3.1) 

j=o-,7r 

where Vq is the fermionic mean field term 

Po(^o(T),T) = - [Vo(^o(T)) - Vo(^o(0))] + ^ j^xln [l + g-/^^^ + V^o(T) ^ (3 3) 

The Vi terms, which contain the contribution of the bosonic fiuctuations, can be divided 
for later convenience into two terms 

Vi ^ + Vi^ (3.3) 

where the subscripts E and M remind that the bosonic propagator is evaluated for Euclidean 
and Minkowski four-momenta respectively. Their explicit forms are 



/■+00 r+00 



p+00 n+00 

/ p dp du \n 
Jo Jo 



{2'KfJo ^ ^io 



\ ^ ^P/3](^^,p;<^o(T),T) 
^^^a)(^^5P;<^o(^)) 



(3.4) 
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with Vi given in Eqs. (^) and 

1 



' IM 



- lim ^ J^P 1^ nB{uj) [arg (iP^. ^{uj + ie, p; v^olT"), T) 



(3.5) 



In the last expression, which is a pure temperature correction, nB{oj) is the Bose-Einstein 
distribution function 

1 



(3.6) 



e/^^ + 1 

and arg{z) is the argument ^ of the complex number z = |2;|exp(i^). The expressions for 
the complete inverse bosonic propagators at finite temperature are given in Appendix 

E 

Notice that for a correct expansion of the pressure, both Vq and Vi are evaluated 
at the mean field solution <^o(^) of the gap equation and also that the bag constant 



i?^Vo(^o(0),0) + - V^(^o(0),0) 

j=(T,7r 



(3.7) 



has been subtracted from the beginning in order that Po = "^i = ^ = at T = (see Eqs. 
(P) and (U)). 



Let us first recall the results |ll^ , p!0[| which are obtained by extremizing the finite tem- 
perature effective potential in the mean field approximation, whose expression is given by 
minus the pressure in Eq. ( p.2|) and with (poiT) taken, in this free variable. 

The gap equation, in D = 2 + 1, can be solved analytically. One obtains 



gVo(y.,T) _ 
dip vr 



m 



mo 



+ a 







(3.8) 



Thus it is immediate to find that for a = 0, the zero temperature solution outside the 
origin <^o(0) = mQ merges continuously into the one located in the origin at the second order 
critical temperature 

mo 



21n2 



0.7213 mo 



(3.9) 



Of course, the correction spoils the validity of chiral symmetry breaking in the case 
of a continuous symmetry ( |2.3| ). As already said, we have considered the latter only for 
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In Eq. (|3.4| ) we have separated the zero temperature effective potential (which depends 
on temperature only through (^o(^)) from the second term, which depends explicitly on 
temperature. 

Apparently Vl^ become ill defined above a certain value of the temperature, since f^oiT) 
decreases for increasing temperatures and, as already said, both and in ( f^.8D and ( f^.9[ ) 
exhibit an imaginary part for small arguments. Furthermore, also the finite temperature 
corrections to them (the logarithms in Eq. ( p.4D ) turn out to be complex. Nevertheless we 
have verified that their sums, which can be written as follows 

He = "■^^^ X P X [^^o^i (^^' P; MT)) + iT^pjii^, P, <^o(^), + c.t. (3.10) 

are always real when evaluated at the minimum (po{T) of the mean field effective potential. 
This is a very important feature since it ensures that the total pressure is a well defined 
quantity even at order 

Finally, the V{j^ terms in Eq. ( p.5|) are pure temperature terms, whose interpretation is 
more direct, at least in case of a free propagation, where they correspond to the pressure of 
a gas of free bosons. 



IV. NUMERICAL RESULTS 
A. Discrete Symmetry 

By handling the expressions for tD^^^ as shown in Appendix |^, the number of successive 
integrations necessary to calculate the pressure in Eq. ( pTT|) is reduced, and thus the study 
can be better completed numerically by means of suitable routines of integration. The 
preliminary treatment of is important expecially to calculate V{^^ which are slowly 
convergent expressions for large momenta. 

Figs. summarize the behaviour of the total pressure for the discrete symmetry model 
(As = 0) in the chiral limit M = (a = 0). With only the scalar composite included, the 
total pressure in Eq. (|3Tl| ) is composed of three terms 
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V = Vo + j^VI = Vo + - {V^E + n*/) (4.1) 

The results are given for A = vr, thus mo = 1 and from Eq. (p.9|) the mean field critical 
temperature is Tc = (2 In 2)~^. Since Tc sets the natural scale of low and high temperatures, 
all the temperature behaviours are plotted vs. T/Tc. As usual, in Figs. the pressure 
is divided by = in order to make evident the Stefan-Boltzmann regime at high 
temperatures. In Fig. |l] the mean field term Vq (stars) is compared with V^^ (empty circles) 
and Pfi^ (full circles) separately. Notice that while V^j^ remains smaller than the mean field 
term in the whole range of temperatures, V^^ is large for intermediate temperatures (even 
if it has still to be divided by A^), and it is slightly negative for T ^ 1.7 T^. In Fig. ^ the 
mean field term is compared to the total contribution of order 1/A^ for = 1, "Pf, which 
comes out to be positive in the whole range of temperatures. 

In Fig. 1^ we plot the total pressure of Eq. ( [4.1|) divided by T^, for A^ = 3 (full triangles), 
A^ = 10 (full circles) and N = oo (stars). The curves show that in absence of pion-like 
excitations, the low-temperature behaviour is not much affected by the 1/A^ correction, as it 
was expected, whereas in the region of intermediate temperatures (around T^) this correction 
is still large with A^ = 10, suggesting the importance of higher orders in 1/A^. Finally, we 
notice that the contribution of the scalar composite starts to be negligible at about twice 
the mean field critical temperature. 

B. Continuous Symmetry 

In Figs. we plot the behaviour of the pressure divided by (see Eqs. ( p.l|) and (|3l3|)) 
for the continuous symmetry model (A5 = A), for nonvanishing fermion masses {a = 0.001 
in Figs. H-^a) and a = 0.01 in Figs. ^-§|(b)). All the curves are plotted vs. T/T^ as done in 
the discrete symmetry case [T^ = (2 In 2)^^). 

In Fig. ^ we compare the mean field term Vq (stars) with (empty circles), Vij^j (full 
circles), PJg (empty triangles) and V^p^ (full triangles). This figure show that, for very low 
temperatures and small values of a, the pressure is dominated by the pion contribution. We 
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will further comment on this point when referring to Fig. |^. Notice also that Vlj^.j and Vfj^ 
become almost degenerate for T > T^. In fact the inverse propagators i'D'^^ and iT)^^ in Eq. 
( p.5[ ) are equal for i^o{T) = 0, and even if, after inserting a bare fermion mass, this condition 
is never realized, for small a, <^o(^) ~ and therefore the difference between Vi^f and Pf^ 
is negligible on the scale of the figures. 

In Fig. ^ the mean field term is compared to the total contribution of and Vi . 

In Fig. 1^ we plot of the total pressure in Eq. for = 3 (full triangles), = 10 

(full circles) and N = oo (stars). 

We remark that, although "Pf is negative in a finite range of temperatures below Tc, the 
total pressure of order 1/A^ is always positive, and that its contribution is still large in the 
critical region, even for = 10. As in the discrete symmetry version of the model the role 
of composites becomes negligible at T ~ 2Tc 

As we have just pointed out. Figs. ||-| show that, for very low temperatures, the pressure 
is dominated by the pion contribution (at least for small values of a). This is in agreement 
with what one expects on general grounds since the contribution of massive particles to the 
partition function is exponentially depressed as exp(— m/T) for T 0, and therefore the 
pion, being the lightest particle in the spectrum, dominates the other terms 

In Fig. I it is also evident that the dominant pion contribution comes in particular 
from the term V^jy.^ given in Eq. (|3.5| ). It is easy to verify that this term reduces to the 
standard expression for the pressure of a free gas, in case of iV^^ = tu^ — — m^. Thus it 
is interesting to compare our results with the free gas approximation in order to quantify, 
in the low temperatures region, the contribution of the pion and sigma poles to the total 
pressure. In the case of free particles the expression for the pressure density is 




(4.2) 



where g is a degeneration factor and ± refer to fermions and bosons respectively. 
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In Fig. 1^ we compare, for various values of a, the low temperature behaviour of Pq, 'Pim 
and Vij^j (divided by T^) with that of a free particles gas of mass m = |(^o(0)|, rn„ and 
respectively (all the masses are evaluated at T = 0). For the fermionic (mean field) term we 
have g = 4: (see Eq. ( |3.2| )) and m = |(^o(0)|, whereas for the bosons the degeneration factor 
is g = 1 and m„ and are defined through the zeros of the inverse propagator in Eq. ( [A3[ ) . 
The curves are plotted for a = 0,0.001,0.01. The pion pressure for a = is plotted as a 
reference level, even though in this model (with continuous symmetry) the zero mass limit 
is not allowed since the pion cannot exist as Goldstone particle in d = 2 spatial dimensions 
at finite temperature. Notice that, whereas the pion term changes dramatically when a is 
turned on, the fermionic and sigma terms for a = 0.001 and a = 0.01 are, on this scale, 
almost undistinguishable from the case a = 0. 

From Fig. ^ it is evident that, for very low temperatures (T ^ 0.2 T^), the data of our 
model are very well fitted by the free gas approximation. This implies that in this region, 
although the boson propagator has a complex structure (see Appendix |^), the pole gives 
almost the whole contribution. For the fermions the approximation works even better due 
to the fact that the only deviation from a free gas behaviour relies on the temperature 
dependence of the condensate (fo, which is almost constant for low T. 

Notice also that, when the exact curves of Vlj^.j start deviating from the free gas be- 
haviour, the contribution of V^^ and is already important (see Figs. §(a)-(b)) and the 
approximation is no longer reliable. In this model this happens for temperatures of the order 
of 0.1 T,. 

V. CONCLUSIONS 

The paradigm of chiral restoration beyond a critical temperature and the ensuing physical 
implications is a crucial aspect of chiral theories, engendering at this time intensive studies. 
The most relevant physical case in high energy is QCD, but its complexity makes its complete 
study still very difficult. For certain aspects, expecially concerning chirality, the Gross-Neveu 
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model is believed to have features similar to QCD. 

We have here performed a calculation of the thermodynamics of the Gross-Neveu model 
in 2+1 dimensions (both in its discrete symmetry and continuous symmetry versions), at 
order 1/N , for the whole range of temperatures below and after the chiral phase transition. 
The calculation is performed by specializing a general formula valid for four-fermion theories. 
Renormalization is treated consistently. 

The possibility of a detailed calculation for this model allows for quantitative comparisons 
of the different contributions to the equation of state in the various temperature regions, 
below, near, and after the transition. In particular we find that the intuitive expectations on 
dominance of the pion-like excitations at low T, and on free-component behaviour at large T 
are satisfied. On the other hand we are able to quantitatively assess the limits for the validity 
of such behaviours. Pion-like dominance holds only for T much smaller than the critical 
temperature, and free-component behaviour only from at least twice the critical temperature. 
Such hmitations hold even for sufficiently large N. Within the intermediate region around 
the phase transition temperature, the role of composite states appears important, and near 
the critical temperature it might be unavoidable to go to higher 1/A^ orders. We also find 
quantitative differences between the cases of vanishing and non- vanishing fermion masses. 
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APPENDIX A: BOSONIC PROPAGATOR 

In this Appendix we discuss the structure of the inverse bosonic propagator, which can 

be written as the sum of the zero temperature term, iV^^, and a temperature dependent 
part, iVp^ (for convenience we omit the index j) 
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(Al) 



The zero temperature inverse bosonic propagator Wq ^ is |T^J§] 



iVo\uj, p; ^) = -(1 + 5Z2) + AiX J 



(A2) 



(27r)3 [(g^ - <^2] (^^^^ - y,2) 

where ± refer to the scalar and pseudoscalar fields respectively, = {u, p) and SZ2 is the 
counterterm needed for the renormalization of the one loop effective potential (see Sect. y). 
The renormalized expression is 



A 



Wq ^ {to, p, ip) = -6{p'^ - uj'^)<mo - \ip\ + 

TT 



1 - p2 



'M 
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tan" 



-1 p'^-uj^ 



(A3) 



A 



+ -e{uj^ -p'^)\mo~\ip\ , ^ 

TT I 4 VCU^ — p^ 



1 / J2 ^2 

iuj - p - ej,,j 
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<^ 


+ Vcu^ - p^ 


2 




- v/t^2 _ p2 



2^ TT ^2 _ p2 _ 



TT ^ ^4 ^yUJ^ — p^ 



where e^^ = 4<p'^ and e^^^ = 0. 

At finite temperature, the inverse propagator given in Ref. [^, can be cast in a more 
compact form as 



iV^\u,p;ip,T)=A\ J 



(fq UFjEq) 
(27r)2 Eq 
(Pq np{Eg 



(27r)2 E, 



A (u^ - p' 



E 



9 a,e2=± 



, (E, + eia;)^ - 



(A4) 



where Eg = vqM-^, and np[Eq) is the Fermi-Dirac distribution function 



(A5) 



To calculate the real and imaginary part of Eq. ( [A^ ) we evaluate iV^ ^ (cu + ie, p; yj, T) in the 
limit e ^ 0. Therefore, by using the formula 

1 



hm / d9- 



e^o Jo A + B COS 9 + iae 

\sign(A)e ( - B^) - isign(a)e (b^ - A' 



(A6) 



we get (x = Eg) 
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with 



,2 „2 ,2 \ 



^ -P 



\A^-B^\ 



sign{A)9 (A^ - B^) - isign(c<; + ^ix)9 [B^ - A^) \ (A7) 



A= (uj'^ - p2) + 2^iu;x 



(A8) 
(A9) 



By studying the sign of — B"^ in Eq. ( |A9| ) , the 6'-functions can be implemented in the 
integration extrema as follows 



/ dx9 (a^ -B')=9 (p2 -cu^) ' dx + 9 ({cu^ - p'){p^ - cu^ + V)) / dx 



2 



+^ (o;^ - p^ - 4y.^ 



+ 00 



(AlO) 



+P\/A 



/^l - = ^ (P' - U-,J- + 9{u^-p'- 4^^) 5,,,., (All) 



where we have defined A = 1 — 4(/)^/ (cj^ — p^) 

The expression of the inverse propagator evaluated for imaginary values of the energy, 
uj — > iuj, needed in Eqs. ( p.4|) and (|3.5| ), follows from Eqs. ( |A3|) and ( |A4|) . At zero tempera- 
ture, from Eq. (|A3|), we easily get 



A 



W^^ {iu,p-^) = -\mo- ^ - - — J 

TT I 2 VCU"^ + p^ 



4(^2 



(A12) 



whereas at finite temperature, by performing the angular integration in Eq. (^), we have 

+ — 6^ — 



iVa^ {iiu,p; (p,T) = — dx np{x 



TT J|(p 



72 (cj2 + p2 + e 



-2 



(A13) 



with 



= (u;^ + p 
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¥ = 4p^q^ = 4p^x^ - 4pV 
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(AM) 
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FIGURES 

FIG. 1. Plot of the pressure entering Eq. divided by vs. T/T^ {T^ = (21n2)"i) for 

the discrete symmetry model (A5 = 0) with zero current fermion masses. The mean field term 
Vq /T^ (stars) is compared to the correction for = 1, which is divided into two terms: 
'Pf^/T^(open circles) and 'Pfj^^/T'^ (full circles). The low temperature behaviours are enhanced in 
the right upper part of the figure. 

FIG. 2. Plot of the pressure divided by vs. T/T^ in Eq. (|1]) (T^ = (21n2)-i), for the 
discrete symmetry model with zero current fermion masses. The mean field "Pq/T^ (stars) is 
compared to the correction for = 1, Vi /T^ (full circles). The low temperature behaviours 
are enhanced in the right upper part of the figure. 

FIG. 3. Plot of the total pressure divided by vs. T/Tc in Eq. (^) (T^ = (21n2)-i), for 
the discrete symmetry model with a = 0. The curves represent V/T^ for = 3 (full triangles), 

= 10 (full circles) and N = 00 (stars). The low temperature behaviours are enhanced in the 
right lower part of the figure. 

FIG. 4. Plots of the pressure divided by vs. T/T^ in Eqs. (|]|), (p|) (T^ = (21n2)"i), 
for the continuous symmetry model with a bare fermion mass term. Figs. (a,b) are for a = 0.001 
and Q = 0.01 respectively. The curves represent Vf^/T^ (open circles), Vf^j/T^ (full circles), 
Vo/T^ (stars), V^f^jT-' (fuh triangles) and Pf^j/T^ (open triangles). V{m/T^ and V^m/T^ be- 
come (almost) degenerate after Tc, the difference being undistinguishable on this scale. The low 
temperature behaviours are enhanced in the right upper part of the figure. 

FIG. 5. Plots of the pressure divided by vs. T/Tc in Eq. {Tc = (21n2)~i), for the 

continuous symmetry model with a bare fermion mass term. Figs. (a,b) are for a = 0.001 and 
a = 0.01 respectively. The curves compare the mean field approximation Vq/T^ (stars) with the 
total 1/A^ corrections for A^ = 1 of the scalar composite Vi /T^ (full circles), and of the pseudoscalar 
VI I T'^ (full triangles) The low temperature behaviours are enhanced in the right upper part of the 
figure. 
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FIG. 6. Plot of the total pressure divided by vs. T/T^ in Eq. (^) (T^ = (21n2)-i), for 
the continuous symmetry model with a bare fermion mass term. Figs. (a,b) are for a = 0.001 and 
a = 0.01 respectively. The curves represent V/T^ for = 3 (full triangles), iV = 10 (full circles) 
and N = oo (stars). The low temperature behaviours are enhanced in the right lower part of the 
figure. 

FIG. 7. Comparison between the low temperature behaviour of Vq/T"^ (stars), Vij^j/T^ (full 
circles) and Vfj^j/T^ (full triangles) and that of a free particles gas of masses m = 
and respectively (dotted lines) for a = 0,0.001,0.01 {T^ = (2 In 2)^^). The pion pressure for 
a = is plotted as a reference level, even though in this model (with continuous symmetry) the 
zero mass limit is not allowed. Notice that, whereas the pion term changes dramatically when a 
is turned on, the fermionic and sigma terms for a = 0.001 and a = 0.01 are, on this scale, almost 
undistinguishable from the case q = 0. 
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